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The dispersive sweep of fast radio bursts (FRBs) has been used to probe the ionized baryon 
content of the intergalactic medium, which is assumed to dominate the total extragalactic 
dispersion. While the host galaxy contributions to dispersion measure (DM) appear to be 
small for most FRBs?, in at least one case there is evidence for an extreme magneto-ionic local 
environment? and a compact persistent radio source Here we report the detection and 
localization of the repeating FRB 20190520B, which is co-located with a compact, persistent 
radio source and associated with a dwarf host galaxy of high specific star formation rate 
at a redshift z = 0.241 + 0.001. The estimated host galaxy DM ~ 903*77, pc cm^, nearly 
an order of magnitude higher than the average of FRB host galaxies?l$| far exceeds the DM 
contribution of the intergalactic medium. Caution is thus warranted in inferring redshifts 
for FRBs without accurate host galaxy identifications. 


FRB 20190520B was detected with the Five-hundred-meter Aperture Spherical radio Tele- 
scope (FAST)! in drift-scan mode as part of the Commensal Radio Astronomy FAST Survey 
(CRAFTS)® at 1.05-1.45 GHz in 2019. Four bursts were detected during the initial 24s scan. 
Monthly follow-up tracking observations between 2020 April and 2020 September detected 75 
bursts in 18.5 hrs with a mean pulse DM of 1204.7 + 4.0 pc cm ?. Assuming a Weibull distribu- 
tion of the burst waiting time, we model the FRB burst rate to be R = 4.5*} 2 hr !, for a fluence 
lower limit of 9 mJy ms and a burst width of 1 ms, indicating that the FRB can episodically have 
a high burst rate. Similar to other repeating FRBs, this FRB shows complex frequency-time in- 
tensity structure?! with multi-component-profiles, sub-burst drifting, and scattering and 
Methods). We detected no linear polarization for FRB 20190520B from FAST observation (see 
Methods), the rotation measure (RM) for this source have been detected from higher frequency 
band“ The properties of FRB 20190520B can be found in Table 1] 


We localized FRB 20190520B with the Karl G. Jansky Very Large Array (VLA) using the 
realfast fast transient detection systent2 Throughout the second half of 2020, we observed the 
source for 16 hours and detected 3, 5, and 1 bursts in bands centered at 1.5 GHz, 3 GHz, and 
5.5 GHz, respectively. We measured a burst source position in the International Celestial Refer- 
ence Frame (ICRF) of (R.A., Decl.) [J2000] = (16h02m04.272s, —11?17'17.32") with positional 
uncertainty (107) of (0.10", 0.08") dominated by systematic effects (see Methods). Deep images 
using data from the same observing campaign revealed a persistent radio continuum counterpart at 
(R.A., Decl.) [J2000] = (16h02m04.261s, —11?17'17.35") with 1c position uncertainty of (0.10", 
0.05") that is compact («0.36"; see Methods) and has a flux density of 202 + 8 wy averaged over 
span of ~ 2 months from August 30 to November 16 of 2020 at 3.0 GHz. Using the average flux 
density of each sub-band over the VLA campaign, we find that the radio continuum counterpart 
spectrum can be fit with a power-law spectral index of —0.41+0.04 (see Methods). 


In order to identify the host galaxy, an optical image (R’-band) was obtained using the 
Canada-France-Hawaii Telescope/MegaCam. [Figure 2]shows the FRB location compared to deep 
optical, near-infrared (NIR), and radio images of the field. The optical image (R’-band) reveals the 
galaxy J160204.31—111718.5 at the location of the FRB, the light profile of which peaks at ~ 1" 
south east. Given the measured offset of 1.3" between the FRB and the peak of the galaxy light 
profile, and the sky surface density of galaxies with this magnitude, we estimate a chance coinci- 
dence probability of 0.8% (see Methods), indicating that J160204.31—111718.5 is the host galaxy 
of FRB 20190520B. The NIR image (J-band, 1.153—1.354 um) obtained with Subaru/MOIRCS!3! 
most likely shows the stellar continuum emission of the host galaxy with an AB magnitude of 
22.1+0.1 in J-band, with FRB 20190520B and the radio continuum counterpart on the galaxy 


periphery. 


We obtained an optical spectrum at the location of the FRB with the Double Spectrograph 
on the Palomar 200-inch Hale Telescope, which reveals the redshift of the putative host to be 
z = 0.241 + 0.001 based on a detection of strong Ha, [O III] 4859À, and [O I] 5007A lines 
(see Methods). A follow-up observation with the Low Resolution Imaging Spectrometer (LRIS) 
at the Keck I Telescope covering both the FRB location and the nearby Subaru J-band source 
along the extended R’-band structure indicates that the R’-band structure is dominated by the [O 
III] emission at the same redshift of z = 0.241. The Ha luminosity Lye = 7.4 + 0.2 x 10 erg 
s! after extinction correction implies a star formation rate of ~ 0.41 Mo yr !. Based on the J- 
band magnitude, we estimate the stellar mass of the host galaxy to be ~ 6 x 10* Mo. Thus, we 
characterize J160204.31—111718.5 as a dwarf galaxy with a relatively high star-formation rate for 


its stellar mass compared with local SDSS galaxies At the luminosity distance of 1218 Mpc 
implied by the redshift, the radio continuum counterpart has a spectral luminosity of L5 Gg; = 
3 x 10? erg s! Hz. 


The redshift of an FRB source is often estimated from the DM attributed to the intergalactic 
medium (IGM), DMigm = DMrggs - DMpost - DMyw. Theoretical calculationsP! and observations! 
have independently estimated the IGM contribution to FRB DM as a function of redshift 
B). For a DM contribution from the Milky Way of 113 + 17 pc cm? (including a +40% range 
for the NE2001 estimate! of the disk contribution and a uniform halo distribution from 25 to 
80 pc cm?) combined with an assumed host galaxy DM of only 50 pc cm ^, the implied redshift 
range for FRB 20190520B is z ~ 2.2 to 0.9 for baryon fractions figm of 0.6 to 1 for the ionized 
IGM, much larger than the measured value (see Methods). 


Instead, using DMigm(z = 0.241) =~ 19572)? pc cm^? (68% interval to account for cos- 
mic variance of DMjgm and 0.85 for the ionized baryon fraction) combined with the Milky 
Way contribution, a very large host DM (disk + circumgalaxy + local to source contribution) 
of DMhost = 903*77, pc cm”? is inferred from the posterior distribution. Over a broader range of 


ionized fractions from 0.6 to 1, DMpost extends from 1020 to 745 pc cm ^. 


In addition to the low chance coincidence probability, the measured DM and scattering prop- 
erties of FRB 20190520B render it unlikely that J160204.31—111718.5 is a foreground galaxy 
with the true FRB host much further in the background. The observed Ha emission implies a DM 
contribution from warm gas in the galaxy ranging from 230 to 650 pc cm? (observer frame) for 
temperatures from 0.5 x 10^ to 5x 10^ K (see Methods). The scattering contribution of a foreground 
galaxy depends not only on its DM contribution but also on the geometric leverage effect, which 
will increase the scattering by several orders of magnitude relative to scattering in the host galaxy. 
If J160204.31—111718.5 were a foreground galaxy, we estimate a fiducial range of scattering times 
from 0.6 to 20 s at 1.25 GHz, orders of magnitude larger than the observed mean scattering time 
of 10 +2 ms. While this estimate depends on parameters such as the ionized gas temperature, path 
length through the galaxy, and the degree of turbulence in the gas, which are not well constrained, 
the observed scattering time is likely too small for the proposed host galaxy to lie in the foreground 
(see Methods). 


The large DMpost inferred for FRB 20190520B demonstrates that the distribution of DMhost 
values for the FRB population has a long tail, which may add considerable variance to estimates 


for the IGM contribution. It is conceivable that the DMpost distribution may differ for repeating 
and non-repeating FRBs, which could make non-repeating FRB DMs more accurate proxies for 
redshift. To understand this further, a larger sample of precisely localized FRBs with measured host 
galaxy redshifts is needed to statistically characterize host galaxies, their circumgalactic media, 
and near-source environments along with the IGM” Accordingly, searches need to accommodate 
large values of DM) os as part of the DM budget for FRB sources. Large DMpost contributions may 


also imply large scattering, which can also reduce search sensitivity. 


The co-located radio continuum counterpart, the star-forming dwarf host galaxy, and the high 
repetition rate make FRB 20190520B a clear analog to FRB 20121102A, the first known repeat- 
ing FRB“land the first to be identified with a compact, luminous persistent radio source (PRS)! 
(Lis cuz & 3x10? erg s^! Hz !). Another repeating source, FRB 20201124A, was also associated 
with a radio continuum counterpart 29:22? however, optical spectroscopy and radio interferometric 
measurements demonstrated that the persistent radio emission was spatially extended and consis- 
tent with an origin in star formation in the host galaxy”! By contrast, the continuum counterpart 
to FRB 20190520B appears to not be from star formation because its luminosity would imply a 


star-formation rate of ~10 Mo yr^! 


, a factor of 25 larger than that measured for the host galaxy 
and a factor of five larger than the highest observed star formation rate for galaxies of this mass [14] 
Given its extreme luminosity, unresolved structure in VLA observations, and offset from the peak 
optical emission in the host galaxy, we conclude that the radio continuum counterpart is a compact 
source (« 1.4 kpc ) physically connected to FRB 20190520B, a PRS like that associated with FRB 


20121102A (which has a VLBI-confirmed size < 0.7 p¢!2). 


More than a dozen FRBs were localized prior to FRB 20190520B, including five repeating 
sources 4P5| but only FRB 201211024 had been associated with a compact PRS. FRB 20121102A 
also demonstrates a sporadically large burst rate24 (e.g. a peak burst rate of 122 hr^!), a substan- 
tial rotation measure that varies over both short (burst to burst) and long (year) timescales?H| and 


DMhost as large as 4300 pc cm? 


, suggesting that burst activity may be correlated with both rel- 
ativistic plasma emitting synchrotron radiation and the presence of thermal plasma in the local 
FRB environment. In addition to a PRS, FRB 20190520B shows a DMmy, that is almost three 
times larger than that of FRB 20121102A, and may have a comparably large RM (Methods). If 
a sizeable fraction of DMyos is from thermal plasma in the circum-source medium, then perhaps 
the presence of a PRS and a dynamic magneto-ionic environment are correlated with FRB forma- 
tion, and repeating FRBs like FRB 20121102A and FRB 20190520B are younger sources that still 


reside in their complex natal environments. 


However, until further VLBI observation constraints become available, it is unclear how 
much of the large DMpost for FRB 20190520B is attributable to ionized gas in the circum-source 
medium vs. the host galaxy’s interstellar medium. Moreover, other repeating FRBs have very deep 
limits on PRS counterparts2228| which complicates the connection between burst and magneto- 
ionic activity and PRS luminosity. For example, FRB 20200120E has an upper limit on persistent 
radio luminosity Lı 5 Gu, < 3.1 x 102 erg s^! Hz! 28 It is associated with an old stellar population 
and has a modest DMpos x 50 pc cm? where are in sharp contrast with FRB 20121102A and 
FRB 20190520B. One possibility is that there are multiple kinds of sources that can emit FRBs, a 
point that has been argued based on burst rate and phenomenology29P0| Alternatively, the source 
properties may evolve as the source ages, the PRS fades, the event rate drops, and the surrounding 
plasma dissipate The similarity of FRB 20190520B to FRB 20121102A suggests a potential 
connection between burst activity, the presence of a PRS, and DMy,, for at least some FRBs. 


Various methods have been used to argue that repeating and non-repeating FRBs comprise 
different subclasses?PP?| either at different evolutionary phases or entirely different physical sce- 
narios. While the observed burst repetition and morphology can be time-dependent due to various 
mechanisms22l PRS emission and DMy,, may reflect more persistent aspects of the FRB environ- 
ment and thus may be more reliable tracers of any putative subclasses. Further progress will result 
from better characterization of the full range of host galaxy DMs, which will also help mitigate 
biases in the DM-redshift relation! 
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Figure 1: Bursts from FRB 20190520B, shown as dynamic spectra and burst profiles. Top panel: Frequency 
integrated burst normalized intensities that were detected during FAST observations. These six bursts are chosen from 
different observation epochs. The burst labels *PXX" are in order of arrival time and the corresponding burst properties 
are shown in supplementary information Table 1. Bottom panel: De-dispersed dynamic spectra of the bursts clearly 
showing their band-limited nature. The DMs are taken from the supplementary information Table 1. The color map is 
linearly scaled. and darker patches represent higher intensities. The bad frequency channels are masked and labeled 
using red patches on the left. The time and frequency resolutions for the plots are downsampled to 0.786 ms and 3.91 
MHz, respectively. 


Figure 2: Optical, infrared, and radio images of the field of FRB 20190520B. In each case, the box is 40" in 
size and the 2” crosshairs indicate the best FRB position at R.A. = 16h02m04.272s, Decl. = —11° 1717.32” (J2000). 
North is up and East is to the left. (a) An optical R’-band image obtained by CFHT MegaCam covers 5427A-7041A, 
including redshifted H8 4861À, [OIIIJ4959À, and [OIIIJ5007 À emission lines from the host galaxy. The bright streak 
at the lower left is an artifact caused by a bright star outside the field of view. The inset shows the host galaxy in a 5" 
region centered on the FRB position, as indicated by the yellow circle (0.1” radius, corresponding to the 1c position 
uncertainty). (b) The infrared J-band image by Subaru/MOIRCS shows emission only at the location of the peak of 
the optical light profile of the host galaxy. The inset is a 5" region matching the inset in panel (a). (c) The radio VLA 
image (2-4 GHz) shows a compact persistent source at the FRB location. The synthesized beam is shown as an ellipse 
of size (0792 x 0747) in the left corner. 
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Figure 3: Galactic and extragalactic contributions to the DM observed for 15 FRBs. Theses 15 FRBs are 
with firm host galaxy associations and redshifts. Sources localized on initial detection (*one-off" bursts) and repeating 
FRBs are identified separately. Galactic disk contributions are estimated from NE2001l8 with +20% uncertainties, 
along with an additional halo contribution of 25-80 pc cm^? (full range). Red error bars on each extragalactic DM 
estimate represent a conservative full range uncertainty derived by adding the disk and halo ranges. The expected 
median DM contribution of the intergalactic medium and the inner 68% confidence interval are shown by the orange 
line and shaded region (see Methods). The host galaxy contributions DMpos shift observed values to the right of 
the band of extragalactic DM predicted for the intergalactic medium alone. FRB 20190520B is a clear outlier from 
the general trend, with an unprecedented DM contribution from its host galaxy. FRB 201211024 is identified for 


comparison. 


Table 1: Properties of FRB 20190520B 


Burst Parameters 

Right Ascension (J2000) 
Declination (J2000) 
Galactic Coordinates (J, b) 


16^02704.2725 (0.1”) 
—11?17'17.32" (0.08") 
359°.67, 29°.91 


Number of detections* 88 

Mean total dispersion measure ( pc cm?) 1204.7 + 4.0 
Mean burst width (ms) 13.5x1.2^ 
Mean scattering timescale (ms) at 1.25 GHz 10 + po 
Measured fluence (Jy ms) 0.03 to 0.33 
DMuwyw.aisc DMww halo ( pc cm^?)* 60 , 50 
DMhost ( pc cm?) BST 
Luminosity distance (Mpc) 1218 
Isotropic equivalent energy (10?" erg) 3.6 to 40 


Persistent Radio Source 
Right Ascension (J2000) 


16^02704.261* (0.1”) 


Declination (J2000) —11°17 17.35” (0.05”) 
Flux density at 3.0 GHz (wy) 202 + 84 
Luminosity at 3.0 GHz (erg s^! Hz!) 3 x 10°° 


Size (3 GHz) 


< 0.36” (< 1.4 kpc) 


Host galaxy 

Redshift (z) 0.241 + 0.001 
Stellar mass (Mo)? ~ 6x 108 

Ha luminosity (erg s!) 7.4 + 0.2 10" 
SFR (Mo yr») ~ 0.41 


à Including the FAST and VLA observations. 


c 


o0 fm o 


^ 


These are the mean values from FAST bursts fitted with a Gaussian pulse convolved 
with a one-sided exponential (see Methods). 

The MW electron density model from NE2001 and YMW16. 

Averaged over span of ~ 2 months from August 30 to November 16 of 2020. 

Based on the J-band magnitude. 

Based on the Ha luminosity. 
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Main-text Figure Legends 


Figure 1 Bursts from FRB 20190520B, shown as dynamic spectra and burst profiles. 
Top panel: Frequency integrated burst normalized intensities that were detected during FAST ob- 
servations. These six bursts are chosen from different observation epochs. The burst labels ‘PXX’ 
are in order of arrival time and the corresponding burst properties are shown in supplementary 
information Table 1. Bottom panel: De-dispersed dynamic spectra of the bursts clearly showing 
their band-limited nature. The DMs are taken from the supplementary information Table 1. The 
color map is linearly scaled. and darker patches represent higher intensities. The bad frequency 
channels are masked and labeled using red patches on the left. The time and frequency resolutions 
for the plots are downsampled to 0.786 ms and 3.91 MHz, respectively. 


Figure 2 Optical, infrared, and radio images of the field of FRB 20190520B. In each 
case, the box is 40” in size and the 2” crosshairs indicate the best FRB position at R.A. = 
16h02m04.272s, Decl. = —11^?17'17.32"(J2000). North is up and East is to the left. (a) An 
optical R’-band image obtained by CFHT MegaCam covers 5427A-7041A, including redshifted 
Hg 4861À, [OIIIJ4959À, and [OIIIJ5007À emission lines from the host galaxy. The bright streak 
at the lower left is an artifact caused by a bright star outside the field of view. The inset shows 
the host galaxy in a 5" region centered on the FRB position, as indicated by the yellow circle 
(0.1" radius, corresponding to the 1o position uncertainty). (b) The infrared J-band image by 
Subaru/MOIRCS shows emission only at the location of the peak of the optical light profile of the 
host galaxy. The inset is a 5" region matching the inset in panel (a). (c) The radio VLA image 
(2-4 GHz) shows a compact persistent source at the FRB location. The synthesized beam is shown 
as an ellipse of size (0792 x 0747) in the left corner. 


Figure 3 Galactic and extragalactic contributions to the DM observed for 15 FRBs. 
Theses 15 FRBs are with firm host galaxy associations and redshifts. Sources localized on initial 
detection (*one-off" bursts) and repeating FRBs are identified separately. Galactic disk contribu- 
tions are estimated from NE2001"4 with +20% uncertainties, along with an additional halo contri- 
bution of 25-80 pc cm ^? (full range). Red error bars on each extragalactic DM estimate represent 
a conservative full range uncertainty derived by adding the disk and halo ranges. The expected 
median DM contribution of the intergalactic medium and the inner 6896 confidence interval are 
shown by the orange line and shaded region (see Methods). The host galaxy contributions DMhost 
shift observed values to the right of the band of extragalactic DM predicted for the intergalactic 
medium alone. FRB 20190520B is a clear outlier from the general trend, with an unprecedented 
DM contribution from its host galaxy. FRB 201211024 is identified for comparison. 


Methods 
Observations 


FAST CRAFTS is a multi-purpose drift scan survey conducted with FAST using a 19-beam re- 
ceiver operating at 1.05-1.45 GHz, deployed in May 2018 and conducting blind FRB searches 
using multiple pipelines?! FRB 20190520B was discovered on November 16th 2019 in archived 
CRAFTS data which is in 1-bit filterbank format with 196 us time resolution and 0.122 MHz fre- 
quency resolution. In this first discovery observation, 3 bursts were detected in 10 seconds, and 
another burst was detected 20 seconds later. These 4 bursts from the drift scan survey gave a pre- 
liminary location for the source within a ~ 5 arcmin diameter region. Taking the pointing location 
from FAST, 2 follow-up observations were performed with FAST on April 25th and May 22nd in 
2020 with 19-beam mode in which 15 bursts were detected. A monthly observation campaign was 
then conducted by FAST using the ~ 100 mas localization from VLA (next section). After some 
regular telescope maintenance, 10 observations were performed spanning from July 30th 2020 to 
September 19th 2020 in which 60 more bursts were detected. 


Bursts were detected from FRB 20190520B in each FAST monitoring observation. We list 
the properties of those bursts detected by FAST in supplementary information Table 1. The burst 
arrival time is in MJD format and has been transformed to the arrival time at the solar system 
barycentre (SSB) at 1.5 GHz with the DM values from supplementary information Table 1. The 
observed DM,,; is measured using the method from ref. 9! and code from ref. 241 We use the DM 
value that maximizes the structure of the highest S/N burst on a given date, and report this same DM 
for all other bursts detected on that date. Apparent epoch-to-epoch DM variations may simply be 
the result of a combination of variable effects, including intrinsic pulse structure, time-frequency 
drift, flux distribution in frequency, and scattering (which affects bursts chromatically). 


The FAST data were searched using a Heimdall based pipeliné*3!51 For the FRB blind 
search in the 19-beam drift scan survey, the polarizations were added and only the total intensity 
(Stokes I) was recorded. The trial DM range was from 100 to 5000 pc cm? and we matched the 
pulse width using a boxcar search up to 200 ms. The candidates above a S/N of 7 and present in less 
than 4 adjacent beams were manually examined for further inspection. After we identified FRB 
20190520B as a new FRB source, the follow-up burst search was done with a narrower DM range 
(100-2000 pe cm™). Following the VLA localization, the tracking observations only recorded 
data from the central beam but with all the Stokes parameters and higher time resolution (~ 50 us). 
Candidates with S/N>7 correspond to a fluence threshold of 9 mJy ms for a burst width of 1 ms. 


The pulse widths were estimated by a Gaussian profile fit if the burst showed no evidence of 
scattering (see section on scattering below). A sub pulse is recognized if the bridge between two 
closely spaced peaks drops more than 5c” below the higher peak. If the burst shows a scattering 
tail, the pulse width is derived from a combined fit for the Gaussian width and scattering time. We 
roughly estimate the bandwidth of each burst by dividing the whole bandpass into 50 MHz sub- 
bands and identifying the subbands containing burst emission. The burst fluences were determined 
using the bandwidth, temporal width, and amplitude of each burst. 


VLA Following the FAST detection of bursts from FRB 20190520B, VLA observations were per- 
formed from July-November 2020 (DDT project 20A-557) at the most reliable position determined 
using FAST (good to 5 arcmin). Most of the observations were in the B array configuration (with a 
maximum baseline of 11.1 km), with the exceptions that the array configuration was BnA on MJD 
59161 and BnA->A on MJD 59167 and 59169. In total, 11.4 hrs were spent on the source at mul- 
tiple bands. The total bandwidths at L (1.5 GHz), S (3 GHz), and C (5.5 GHz) bands were 1024, 
2022, and 2022 MHz, respectively, with 1024 channels, corresponding to channel bandwidths of 
approximately 1, 2, and 2 MHz. 


The details of the observations are given in Extended Data (ED) Table 1. The telescopes were 
pointed at the field centered at (R.A., Decl.)[J2000] = (16h02m01s, —11d17m28s). We used the 
realfast search system at VLA to search for bursts from FRB 20190520B in our VLA observations. 
Observations on MJD 59169 contributed to the imaging at S band, but the realfast system was not 


run on this day due to a system error. 


The realfast search system has been described in detail in ref¢!2136, but here we discuss it 
briefly. Using a commensal correlator mode, the visibilites are sampled with 10 ms resolution and 
distributed to the realfast GPU cluster to search for transients. The search pipeline rfpipe then 
applies the online calibration, de-disperses the visibilities, and forms images at a set of different 
temporal widths. The 8c fluence limit of a 10 ms image is 0.29, 0.18 and 0.13 Jyms at L, S 
and C band, respectively. Candidates with image S/N greater than the threshold triggered the 
recording of fast sampled visibilities in windows ranging from 2 to 5 seconds centered on each 
candidate. For each candidate, time-frequency data corresponding to the maximum pixel in the 
image are then extracted and post-processed. The candidate is then classified by fetch, a GPU 


based convolutional neural network’! 


Promising candidates selected visually by the realfast team go through an offline analysis to 
refine the candidate parameters and improve its detection significance, if possible. Several methods 
are tried within this analysis: offline search for the transient using a finer DM grid, varying RFI 


flagging thresholds, changing image size (as the realtime-search is done on non-optimal image 
sizes due to computational limitations), sub-band search etc. (see Section 2.4 of ref. = for more 
details). The pixel sizes of images formed by the realtime system at L, S, and C band were 0.9", 
0.48", and 0.27", respectively, with an image size of 2048x2048 pixels corresponding to field 
sizes of 0.5°, 0.27°, and 0.15°. During refinement, we searched the data with smaller pixel sizes: 
0.5" x0.75" at L band, 0.38" x0.28" at S band, and 0.27"x0.18" at C band. The significance of 
an astrophysical transient improves when the data is de-dispersed at a DM closer to the true DM 
of the candidate (see Section 2 of ref. B9, Therefore, we re-ran the search with a finer DM grid 
at 0.1% fractional sensitivity loss, as compared to 5% used by the realtime system. Noise-like 
events or RFI are sensitive to RFI flagging and image gridding parameters, and so they cannot 
be reproduced on refinement and are discarded. Sometimes the transient signal is present only 
in a fraction of the whole frequency band, so we also re-run our search only using the relevant 
frequencies, which further improves the detection significance. We applied these techniques on all 
the candidates selected from these observations to refine their parameters. 


Localization of bursts 


Calibration The realfast system makes several assumptions during calibration and imaging in 
order to improve computational efficiency for the real-time imaging and search. Moreover, the 
Point Spread Function (PSF) of the interferometer is not deconvolved from the image, so the real- 
time system forms “dirty" images. The PSF shape makes the images more difficult to visually 
interpret and model. 


To address these issues, we used a top-hat function to select the raw, de-dispersed, fast- 
sampled visibilities containing only the burst signal to re-image the burst data using the Common 
Astronomy Software Application package (CASA v5.6.2-3¥ Observations of 3C 286 (acquired 
before the FRB observations) were used to calibrate the flux density scale, bandpass, and delays. 
Complex gain fluctuations over time were calibrated with observations of calibrator J1558—1409. 
We performed phase-reference switching on intervals of 16, 13, and 12 minutes for L, S, and 
C-bands, respectively, consistent with the nominal phase-coherence timescales for the VLA, so 
that after calibration any short-timescale phase variations are negligible. Therefore, the systematic 
errors for the burst positions on short timescales are of the same magnitude as those for the deep 
imaging, as discussed below. 


Determining properties of individual bursts After calibration, the CASA task tclean was used 
to generate an image for each burst and estimate the S/N. Most of the bursts are spectrally confined, 


so for each burst we select the spectral window range that produces the highest image S/N. We then 
use the CASA task imfit around the FRB position to fit an elliptical Gaussian to the source in the 
radio image and measure the centroid location, peak flux density and 1o image-plane uncertainties. 


Hereafter, we refer to these image-plane uncertainties as statistical errors on burst positions. 


The average burst positions are calculated for each frequency band separately, by weighting 
each detection by the inverse of the position fit errors (statistical errors) reported by CASA. Statis- 
tical errors are inversely proportional to S/N, and therefore high significance detection is expected 
to have smaller fit errors. The total positional error at each frequency band is obtained by adding 
this statistical error in quadrature with the systematic error in that band (determined using the deep 
radio images described below). The total positional error (R.A. error, Decl. error) is (0.25", 0.32") 
at L band, (0.28", 0.17") at S band, and (0.12", 0.09") at C band. In the final analysis, the errors are 
dominated by systematic uncertainties at each frequency band (see [ED Table 2|and [ED Table 3). 
The final burst position was obtained by taking a weighted average of the burst positions at each of 
the three frequency bands, using the inverse of the total error at each frequency band as the weight. 
The best estimate of the burst position is R.A.= 16h02m04.272s, Decl.= —11?17'17.32"(J2000). 
We estimate the error on this position to be 0.10" and 0.08". We calculated the reduced chi-square 
value of individual burst positions with respect to this best estimate of the FRB position. The re- 
duced chi-square was given by x? = (1/8) Y; (60; — Bs / Up. where 6; refers to the coordinates of a 
burst, Ovest refers to the best estimate of FRB position, c; is the total error on a burst (statistical and 
systematic errors added in quadrature) and the sum is over all the 9 localized bursts. We obtained 
a reduced chi-square value of 0.53 and 0.77 for R.A. and Decl. respectively. 


We confirmed that there were no short-timescale phase errors by performing intermediate- 
timescale imaging of the continuum data. We imaged segments of 5-30 seconds, such that at least 
1-3 sources could be detected at each frequency, and inspected the position variations of those 
sources over time. At all frequencies, the positions were stable and consistent with the statistical 
uncertainty (radiometer noise). Offsets of the sources were within the range of the systematics 
quoted for the PRS in next section. 


Due to the low time resolution of VLA data (10 ms), we do not perform any sophisticated 
modeling of burst properties. The properties of the VLA bursts, including the flux densities, were 
obtained using CASA and the realfast system, and are given in[ED Table 4| Here, the reported 
widths should be considered as upper limits. The frequency extent of the burst signal was visually 
determined and reported in the last column. We also could not estimate a structure maximizing 
DM for VLA bursts because of the coarse time resolution. The DMs reported in this table are the 


values that maximize the S/N. This also explains the apparent variability seen in the DM values. 


Persistent Radio Source 


Deep radio images and PRS The VLA campaign obtained two epochs at 1.5 GHz and six epochs 
at 3 GHz and 5.5 GHz, which resulted in an on-source time of ~ 3 hours for 1.5 GHz, and ~ 4 
hours for both 3 GHz and 5.5 GHz. 


Along with the realfast output, the VLA visibilities with 3 s or 5 s sampling times were 
saved and analyzed to search for persistent radio emission. This is done in parallel to saving the 
high time resolution (10 ms) data around the burst that was used in the previous section. We used 
the same data reduction, flagging, and calibration approach that was used for burst localization. We 
then performed further flagging on the target and then subsequently imaged its Stokes I data using 
the CASA deconvolution algorithm tclean. To balance sensitivity while reducing sidelobes from 
a nearby bright source, we imaged with a Briggs weighing scheme (robust=0). In addition, self- 
calibrations were performed for all observations to correct considerable artefacts from the close-by 
bright sources in the field. We made use of the CASA task imfit to measure source flux densities 
by fitting an elliptical Gaussian model in the image plane. 


We stacked observations at each central frequency in the uv-plane and then imaged the Stokes 
I intensity, resulting in deep images at 1.5 GHz, 3 GHz, and 5.5 GHz with rms noise of 9.0, 4.5, and 
3.0 uJy/beam, respectively. At 1.5 GHz, 3 GHz, and 5.5 GHz, the upper limits of the deconvolved 
sizes of the PRS are as large as 1.4"x0.89", 0.51"x0.14", and 0.36"x0.1", respectively. Thus a 
conservative upper limit of the size of the PRS from VLA observations is 0.36", which corresponds 
to 1.4 kpc at the angular diameter distance of 809 Mpc. The obtained positions for the PRS are 


shown in ED Table 2| and the positions of the bursts and PRS are shown in |ED Figure 1| The 


systematic offsets on these positions are estimated in the next section. 


Systematic Offsets In order to determine the systematic errors on the coordinates of the PRS 


that we determined from the deep images, we ran the PyBDSF (https: //www.astron.nl/citt/ 
pybdsf/index.html) package to extract radio sources from the deep images. We then cross- 


matched the detected point sources in the deep images with the sources listed in the optical 
PanSTARRS survey DR 18 We identified the radio point sources using the following criteria: 


e The peak intensity (Jy/beam) of a source should be 0.7, 0.5, 0.5 times higher than its inte- 
grated flux (Jy) for the 1.5 GHz, 3 GHz and 5.5 GHz images, respectively. 


e The S/N (peak intensity / local rms noise) of a source should be greater than 5. 


In total, we detected 375, 113, and 43 sources in the 1.5 GHz, 3 GHz, and 5.5 GHz deep 
image, respectively. The selected sources were also visually checked to make sure that they are 
‘point-like’ sources in the deep images. We searched for matching optical sources within radii of 
0.5", 1.0", and 2.0". Going from 1.0” to 2.0", the additional cross-matched sources at each band 
are consistent with chance coincidence. Therefore, we adopt a 1.0" cross-match radius, finding 
136, 31, 9 sources with optical counterparts in the PanSTARRS DRI catalogue at 1.5 GHz, 3 GHz, 
and 5.5 GHz, respectively. These cross-matched sources were used to determine the astrometry 
for FRB 20190520B and its PRS. By subtracting PanSTARRS coordinates from VLA coordinates, 
we estimate the systematic offsets in the 1.5 GHz, 3 GHz and 5.5 GHz positions, as listed in [ED] 
The systematic offsets are consistent with zero mean, and their uncertainties dominate the 
uncertainties of the PRS position. 


Determining the position of the PRS To determine the positions of the PRS in the three fre- 
quency bands, we followed a procedure similar to what was used to determine the burst positions 
(see Section|Localization of bursts). The best estimate for the PRS position is (R.A., Decl.) [2000] 
= (16h02m04.261s, —11?17'17.35"). The error on this position is estimated to be 0.10" and 0.05" 
or R.A. and Decl., respectively. 


Variability and Spectrum of the PRS The flux density of the source measured at each epoch 
is shown in [ED Table 1| Measured flux densities show variations that are mostly consistent with 
measurement errors, but there are ~ 20 variations in the 2.5 GHz subband that, if real, could be 
refractive interstellar scintillation, or intrinsic variations, or both. In order to study the spectrum of 
the PRS, we split each of the observations into two, 0.5 GHz/1 GHz sub-bands. Then we measured 
the average flux density at each of the subbands over the campaign. The multi-band data were fit 
with a power-law model (S, œ v^, where S, is the observed flux density at the frequency v, œ is 
the spectral index), yielding an average spectral index for the PRS of —0.41 +0.04 (see [ED Figure] 


2). 


Chance coincidence association of the PRS In the 1.5 GHz deep image, we detected 8 ‘point- 
like’ sources, including the PRS, within 5 arcminutes of the phase center and with a flux density no 
less than that of the PRS (at 260 uJy) based on our point-source selection criteria described earlier. 
There is an additional bright source with a flux density of a few tens of mJy in the region, but it 


was not classified as a ‘point-like’ source based on our criteria. 


To estimate the chance coincidence probability of the PRS with the bursts, we compared 


the solid angle corresponding to the uncertainty of burst localization and the average solid angle 
occupied by each of the eight *point-like' sources in the FRB field-of-view. The solid angle corre- 
sponding to each of the 8 ‘point-like’ sources is roughly estimated as S,ource= 7(5/60)?/8 Sr. The 
offset between the average position of the nine bursts and the position of the PRS at 5.5 GHz, 
which is best constrained when taking both statistical and systematic errors into account, is about 
0.06 arcseconds. This, along with a statistical error of 0.01 arcseconds and a conservative estimate 
of the systematic error of 0.12 arcseconds, can be used to estimate the offset between the PRS 
and the FRB position. We conservatively estimate this offset to be 0.19 arcseconds. The solid 
angle corresponding to the offset therefore can be estimated as Sog,4—7 x (0.19/60/60)? Sr. The 
ratio between Sofiset and Syource gives the chance of coincident association of the PRS with the FRB 
position to be ~ 3 x 10. 


Galaxy Photometry and Redshift Determination 


The deep R’-band (5427 À- 1041À) images obtained by CFHT/MegaCam were stacked from archival 
observation data taken in 2014-2015 by the CFHT archival pipeline MEGAPIPE, with a total of ~ 
3.6 hours on the field. The level 3 (flux calibrated) images were retrieved for our analysis. 


NIR J-band images of the FRB 20190520B field were taken under a relatively poor seeing 
condition (^ 1.3") through a Subaru target-of-opportunity observation on August 5th 2020. A 
total of 1.4-hour observations were used for the final combined J-band image shown in [Figure 2] 
A J-band source of 22.07x0.14 mag (AB) was detected at ~ 1" south east of the burst location, 
and is possibly stellar emission from the host galaxy. A faint source 2.5" north has 22.87+40.26 
mag in J-band. Neither of these two sources was detected in the K s-band image, with a 5o limit 
of 21.74 mag (1.1 hours). 


An optical spectrum was obtained with the Double Spectrograph (DBSP) on the Palomar 
200-inch telescope on 24th July 2020 using a 1" slit-width. This observation was executed before 
the CFHT/MegaCAM archival data on FRB 20190520B field were found, and only Pan-STARRS 
images were used for observation planning. The slit of DBSP was set to cover the PRS emission 
at R.A. = 16h02m04.27s; Decl. = -11?17'17.5" detected by VLA in L-band on 22nd July 2020. 
No clear optical counterpart was detected in any of the 5 band images of PanSTARRS from DR1 
Hl The slit was guided by the nearby M-star at R.A. = 16h02m04.48s, Decl. = —11°17’19.1” 
as reported in the PanSTARRS DRI catalog, with i = 20.4 mag, and 3.4" due east of the PRS 
coordinates. The slit was set to a position angle of 108.5 degrees, ensuring that both the PRS and 
the M-star fell within the slit. The observations with 2x900 s exposures were carried out under 


photometric sky conditions and sub-arcsecond seeing. The 2D spectrum was generated using 
IRAF, including bias removal, flat-fielding, and reduction of other instrumental effects. The 1D 
spectrum was extracted from a 1.5" window. The standard star BD+28 4211 was used for telluric 
correction and flux calibration. The DBSP 1D spectrum is shown in [ED Figure 3| The flux scale 
of the spectrum does not include the slit loss and registration error of Pan-STARRS coordinates 
of the M-type star. The [OIII] 5007A line and the Ha line are both well detected (> 5o). The 
two emission lines are narrow, each with a FWHM of ~ 10A. The redshift derived from these two 
spectral lines is z = 0.241+0.001. 


A follow-up Keck LRIS spectroscopic observation was carried out on 25th August 2020 
under reasonable weather and seeing conditions (1.1") . The 1.5" slit was set at a position an- 
gle of 160° to the extended optical emission seen in the MegaCam R’-band image around the 
FRB 20190520B location. A total exposure of 3600s was obtained. The emission lines Ho, Hf, 
[OIII]4859À, and [OIIIJ5007À are well detected, indicating that the extended R’-band structure 
has the same redshift of z — 0.241. After the spectrum was corrected for Galactic foreground 
extinction, the line flux ratio between Ha and Hf was used to estimate an extinction of Ay = 0.80 


mag, assuming a case-B line ratio of 2.86, yielding an Ha flux Fu, = (42.0 + 0.5) x 10 erg s^! 


cm ?. The extinction-corrected Ha luminosity is Lye = 7.4 + 0.2 x 10” erg s™!. 

Chance association probability of the host galaxy We use the approach of ref. to estimate the 
chance coincidence probability of the association between the galaxy (J160204.31—111718.5) and 
FRB 20190520B. Assuming a uniform surface distribution of galaxies, the probability of chance 
coincidence follows a Poisson distribution; i.e., P = 1—exp(—n;), where n; is the mean number den- 
sity of galaxies brighter than a specified R’-band magnitude in a circle of given radius, determined 
by the half light radius of the galaxy and the burst localization error region. This number density 
is estimated using the results from ref. H3 From previous discussion, we assume the localization 
error on the FRB position to be ~ 071, and the size of the host galaxy to be 075. The R’-band mag- 
nitude of the possible host is difficult to estimate since it is significantly affected by the emission 
lines. Thus we conservatively estimate the R’-band magnitude to be < 23.3 mag assuming a flat 
spectral energy distribution in vL, between R’-band and J-band, where v and L, are frequency and 
luminosity, respectively. The probability of a chance association between FRB 20190520B and 
J160204.31—111718.5 is estimated to be ~ 0.8%, independent of constraints from the association 
with a rare PRS or the observed FRB scattering. The next nearest galaxy is 6.5" away from the 
FRB and has a chance coincidence probability of > 20%. 


FAST Burst Sample Analysis 


Repetition rate FRB 20190520B was found to be regularly active with two or more bursts de- 
tected in each monitoring observation with FAST. Because of the possible clustering behavior of 
FRB emission, a Weibull distribution was employed for the waiting time 6 (the separation between 
adjacent detected bursts) Eg 
w= E (è) 
AMA 


where k is the shape parameter and 4 is the scale parameter of the distribution. The Weibull 


k-1 " 
pum (1) 


distribution reduces to a Poisson distribution when k = 1. The mean waiting time, E(ó) = AI(1 + 
1/k), implies a burst rate r = 1/E(6). shows the inferred parameter distributions 
obtained with a Markov Chain Monte Carlo (MCMC). We find the burst rate of FRB 20190520B 
is r = 4.5*12 hr! with shape parameters k = 0.377)" when all 79 bursts above 9 mJy ms are used 
(left panel in[ED Figure 4). Waiting times shorter than 1 s include substructure in individual bursts, 


so we also analyze only longer wait times ô > 1 s, yielding r = 5.3*1) hr ! with shape parameters 
k = 0.7670 53 (right panel in|ED Figure 4). 


Short and long timescale periodicity search For a range of trial periods (P) and period deriva- 
tives (P), the times-of-arrival (ToAs) of the 79 FAST bursts were folded to calculate the phase of 
each burst for a given trial period. For the short timescale periodicity search, the ToAs were folded 
at P between 1 ms and 1000 s and P between 10^? and 1 s s™!. For the long timescale periodicity 
search, the trial range of period P was from 2 to 365 d and the range for period derivative (to ac- 
count for spindown) was P from 10^? to 1 d d^!. The longest contiguous phase segment without 
bursts was calculated for each fold trial as a signature for possible periodicity. 


For the short timescale periodicity search, the bursts spread to a phase window larger than 
60% of one period, which indicates no periodicity pattern in the trial range. For the long timescale 
periodicity search, the longest phase segment without bursts showed two maxima ~ 0.6 near 67 d 
and 169 d. These two maxima are reproduced by folding the MJDs of all observation sessions, 
which indicates that the maxima reflect the sampling window rather than any burst periodicity. 
Thus, no long or short period of FRB 20190520B was detected. 


Energy distribution A 1-Kelvin equivalent noise calibration signal was injected before each ob- 
servation session to obtain a high quality flux density and energy calibration measurement for each 
detected burst. The injected noise calibration signal was used to scale the data into temperature 
units, yielding a nearly constant rms radiometer noise to within 6%. We then converted from 


units of Kelvin to Jy using the zenith angle-dependent gain curve, provided by the observatory 


through quasar measurements 51 For most observations, the zenith angle is < 20 degrees, which 
corresponds to a stable gain of ~16 K/Jy. 


Assuming bursts have spectral index of about 0, we calculated their isotropic equivalent burst 
energies, E, following Equation (9) of ref. 
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where F, = S, x Weg is the specific fluence in units of erg cm ?Hz ! or Jy - ms, S, is the peak flux 


density which has been calibrated with the noise level of the baseline, giving the flux measurement 
for each pulse at a central frequency of v. = 1.25 GHz, W, is the equivalent burst duration, and 
the luminosity distance D; = 1218 Mpc corresponds to a redshift z = 0.241 for the source of FRB 
20190520B, using standard cosmological parameters #7 The fluence-width distribution at 1.25 
GHz for FRB 20190520B bursts can be seen in[ED Figure 5| The histogram of burst energies 
exhibits a bump that we fit with a log-normal function. 


DM inventory analysis The observed DM can be separated into four primary components (all in 


the observer’s frame): 
DMos = DMmw + DMmw,halo + DMicm + DMhost (3) 


where DMyy is the contribution from the Milky Way’s interstellar medium, DMhaio is the contri- 
bution from the Milky Way halo, DMhos the contribution from the host galaxy including its halo 
and any gas local to the FRB source, and DMjgy is the contribution from the intergalactic medium. 


We use the NE2001 model! to evaluate DMyw = 60 pc cm^? (compared to 50 pc cm? 
from the YMW16 model) as the mean of a uniform distribution with a generous +40% width to 
conservatively represent the uncertainty in DMyy for the direction of FRB 20190520B. For the 
MW halo contribution, we use a uniform distribution from 25 to 80 pc cm. Together the MW 
disk and halo components yield a total range from 61 to 164 pc cm? with a mean of 113 pc cm ? 


and RMS uncertainty of 17 pc cm ?. 


For the IGM we use the ACDM cosmological model to calculate the mean DM contribution”, 


DMicm(Z) = 3cHoQb fiom ip | XQ + dz " 
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where the free electron number per baryon in the universe is y(z) ~ 7/8, (assumed constant for 
FRB redshifts), the normalized matter density is ,, = 0.315 + 0.007, the dark energy fraction is 


Q4 = 1 — Qn, the baryonic fraction is Q, = A~? x (0.02237 + 0.00015), the fraction of baryons in 
the IGM is ficm, and the Hubble constant is Hyp = h x 100 kms! = 67.36 + 0.54 kms ! Mpc'! 
E Using values for the speed of light c, the gravitational constant G, and the proton mass mp, the 


resulting expression, 


E z 1 +z)dz 
DMiaw(z) ¥ 978 pc cm? fiom f ( ) : (5) 


0 [Q,(1 +z} + Qa]? 
yields DMiaw(0.241) = 248 fiam pc cm ?. 


Using DMjgm(z) as the mean, we calculate the range of values for a given redshift using a log- 
normal distribution with variance a5... (2) = DMicm(z)DM,, where DM, = 50 pc cm ? is chosen 
to provide cosmic variance consistent with published simulations. This gives &'pw4,(0.241) = 
111 Vftcm pe cm ?. The parameters for the log-normal distribution are then OinDMicm = Un[l + 
(Tpvieu(Z)/DMicm(z))7]}'? and fin pu, = INn[DMicm(z)] — o pneu 2. Note that the log-normal 
distribution is asymmetric. Equation |5| and the log-normal distribution are also used to estimate 
the median DMjcy and its inner 68% uncertainty range shown in 


To constrain the host-galaxy DM, we combine the MW and IGM estimates and their uncer- 
tainties with the measured DM averaged over all bursts, which is DMo»s = 1204.7 + 4.0 pe cm ?. 
The cumulative distribution function (CDF) of the log-normal distribution for the IGM contribu- 


tion yields the range DMigm = 1951? pc cm? 


. We then calculate the posterior distribution of 
DMhos after marginalizing over the IGM and Milky Way distributions and using a flat, uninforma- 
tive prior for DMyos. The median and probable ranges are calculated from the corresponding CDF. 
Using ficm = 0.85 for the baryon fraction in the IGMJ, we obtain a median value for DMhos 
and 68% probable interval DMpos: = 903*77, pc cm. When we vary fiw from 0.6 to 1, the total 
range of values for DMijay is ~ 80 to 350 pc cm ? and for DMpos is ~ 1020 to 745 pc cm 2, where 
we have taken the median values resulting from each value of fig and extended the ranges using 


(half of the) corresponding 68% probable ranges for the smallest and largest values of fiGw. 


For comparison we also give redshift estimates for different values of figm if a fixed value 
of DMnost = 50 pc cm? is used (as is often found in the literature) along with the above quoted 
mean value for the MW contribution. For ficm ranging from 0.6 to 1, we obtain a range of redshift 
values from 2.2 to 0.9 (including the 6896 redshift interval derived from the redshift CDF for each 
value of fjgm), much larger than for the redshift of 0.241 for the identified host galaxy. 


Constraints on DMy, from Ho Measurements The DM of the host galaxy is independently 
estimated from its Hæ emission by converting the extinction-corrected Ha flux, Fu, = (42.0x0.5)x 
107 erg s^! cm ?, to an Ha surface density of 224 +3 Rayleighs in the source frame at z = 0.241, 


assuming host galaxy dimensions of 0.5 by 0.5 arcseconds?!l These host galaxy dimensions are 
only a rough estimate based on the size of the Ha emitting region in the Keck image, but because 
the image is seeing limited the assumed dimensions may be even smaller, which would only serve 
to increase the Ha surface density in the following calculations. The Ho flux and estimated surface 
density are very similar to those found for the host galaxy of FRB 20121102A 511 The Ha surface 
density in the source frame S (Ha), is related to the Ha surface density in the observer frame as 
S (Ha); = (1 + z)*S (Ho, and is related to the emission measure (EM) in the source frame by 


(6) 


S (Ha), 
EM, = 2.75 pc cm TA (Ho), « 616 € 7 pccm ^ x m | (Ha) | l 


224+3R 


where we have used the redshift z = 0.241 and T, is the temperature in units of 10^ K. 
The EM is related to the DM in the ionized cloudlet model =! by EM = [Z(1 + &)/ fIDM? LL, 
where Z represents cloud-cloud variations in the mean density, €? represents the variance of density 
fluctuations in a cloud, f is the filling factor, and L is the path length through the gas sampled by 
the FRB. Using this relation, we obtain the corresponding source frame DM, 


1/2 1/2 
f/0.1 | Peer l (7) 


f(1 + €)/2 224+3R 


where we have adopted fiducial values of f = 0.1 and = e? = 1, which are typical for the 


L yn 
DM, ~ 392 + 3 pc cm? x T2 (s k 
pc 


warm ionized medium in the Milky Way. For the path length through the gas we adopt a fiducial 
value of L = 5 kpc, which is based on the apparent size of the Ha emitting region in the host 
galaxy. However, given that the optical image is seeing limited, and that the orientation of the 
galaxy relative to the FRB line-of-sight is not known, this path length could be as large as 10 kpc 
or as small as 0.1 kpc. For a range of L from 0.1 to 10 kpc, we find that DM, could be as small 


3 or as large as 560 pc cm™?, for the same fiducial values of T4, f, Z, and €?. The 


as 55 pe cm™ 
estimated DM is also highly sensitive to the temperature: for a range of T4 from 0.5 to 5, DM, 
could range from 290 to 810 pc cm ?, keeping all other parameters fixed at their fiducial values. 
In the observer's frame, the measured DM contribution of the host galaxy is smaller by a factor 
of 1/(1 + z), yielding a nominal value of DMhost.coc ~ 316 +2 pc cm ? for the coefficient in 
Equation [7] The quoted errors only account for measurement errors in the Ha flux. To match the 
inferred value of DMpos ~ 900 pc cm? requires that the three factors in Equation [7| involving 
T4, L and f combine to a factor ~ 3, which may easily be explained by a higher temperature or 
the broad range of possible values for L. Regardless, the large Ho EM affirms that the FRB DM 
receives a significant contribution from ionized gas in the host galaxy, but it is unclear whether the 


diffuse, Ha emitting gas can account for the entire host galaxy contribution or whether the FRB’s 


local environment contains significant ionized gas content that is not seen in Ha. 


Extragalactic scattering Scatter broadening manifests as a frequency-dependent temporal asym- 
metry in the burst dynamic spectrum, which is typically modeled as a one-sided exponential pulse 
broadening function convolved with a Gaussian pulse. Some bursts from FRB 20190520B have 
burst structure that is suggestive of scattering: leading edges that are aligned across the radio 
frequency band coupled with temporal asymmetries that broaden at lower observing frequencies. 
However, a number of bursts are symmetric across the radio frequency band and vary significantly 
in burst width. When the burst S/N is low, it is also difficult to distinguish scattering from other 
burst substructure, such as the time-frequency drift of intensity islands?! (often referred to as the 
“sad trombone" effect). Examples of bursts with and without evidence of pulse broadening are 


shown in|ED Figure 7 


To measure the mean burst width and scattering timescale of FRB 20190520B, we first inte- 
grate the dynamic spectrum of each burst along the frequency axis and normalize the resulting burst 
profile. Each burst profile is fit with a model comprised of a Gaussian component convolved with a 
one-sided exponential. Out of the 79 fitted scattering times, all scattering times with fractional un- 
certainties > 50% are excised, leaving a subset of 26 bursts. The scattering times of these 26 bursts 
are then compared to their dynamic spectra and 1D burst profiles to verify that the burst temporal 
asymmetry is broadly consistent with a v ^ frequency scaling. Table 1 in the Supplementary In- 
formation shows the scattering times re-scaled to 1.25 GHz assuming a v * frequency dependence. 
The remaining 53 bursts, whose scattering times have fractional uncertainties > 50%, are then 
re-fit with a Gaussian-only model for the 1D burst profile. The full-width-at-half-maximum burst 
widths are also shown in supplementary information Table 1. This provisional approach yields a 
mean scattering time and burst width of 10 + 2.0 ms and 13.5 + 1.2 ms at 1.25 GHz, respectively. 
A more detailed study investigating the burst widths in separate frequency subbands is in progress, 
and corroborates the scattering interpretation presented here. 


The observed mean scattering time of 10 + 2 ms at 1.25 GHz is likely too small for the FRB’s 
host galaxy to lie behind our proposed host galaxy association. In this alternative scenario, the host 
galaxy would lie at a redshift zp > 0.241 that depends on the foreground galaxy’s contribution to 
the total DM budget. The FRB would pass through the putative intervening galaxy at a redshift z = 
0.241 at an impact parameter of about 4 kpc (based on the observed offset of the FRB localization 
in the optical images). The scattering contribution of the intervening galaxy lens is related to its 


DM contribution by=2/53) 
F C scatt DM; 


(1 t ze)’ ví , 


where DM, is the contribution of the intervening galaxy in pc cm ? in its rest frame, v is the 


t(DM, v, z) ~ 48.03 us x (8) 


observing frequency in GHz, z is the intervening galaxy redshift, and F = £e€?/ f(21;)!? quantifies 


-1/3 


the electron density fluctuations in the scattering layer for F in units of pc? ? km! ^^, where Z, €, 


and f describe the density fluctuation statistics and filling factor, /, is the outer scale of turbulence 
and l; is the inner scale. The fluctuation parameter F in the Milky Way varies from ~ 107? pc-?? 
km-!^ in the thick disk to = 10? pc?? km! near the inner Galaxy, and is typically ~ 0.1 — 1 


-1/3 in the thin disk, but it is not generally known for other galaxies 521 The geometric 


pc ?? km 
factor Gscat is unity for scattering in a host galaxy but it can be very large for an intervening 
galaxy, in which case Gscat & 2d4di,/d,4L, where dy, dio, and dso are the angular diameter distances 
between the source and lens, lens and observer, and source and observer, respectively, and L is the 
path length through the lens. The numerical pre-factor in Equation |8| is for all angular diameter 


distances in Gpc and L in Mpc. 


For the nominal DM contribution implied by the Ha emission, DM; ~ 392 pc cm ?, and a 
path length through the intervening galaxy L ~ 5 kpc, the implied scattering time is tT ~ 20 s at 
1.25 GHz, orders of magnitude larger than the observed mean scattering time. For a smaller gas 
temperature T4 ~ 0.5, yielding a DM contribution DM; ~ 290 pc cm^?, the implied scattering 
time is still r ~ 18 s at 1.25 GHz, assuming F ~ 0.1 pc?? km"! and L ~ 5 kpc. However, 
there is significant latitude in this estimate, depending on the assumed gas temperature, path length 
sampled by the FRB, and the value of F. While we have assumed that the FRB traces the entire 
DM, implied by the Ha emission, the FRB LOS likely only traces a fraction of the Ha emitting 
gas. If DM; is as small as 50 pc cm^, then 7 could be as small as 0.6 s for F ~ 0.1 pc?? km~! 
and L « 5 kpc. While the implied scattering time could be two orders of magnitude smaller or 
larger depending on the combination of Ë, the path length L, and the parameters used to estimate 
the Ha DM, our fiducial estimates suggest that the observed scattering time is significantly smaller 
than the scattering expected from an intervening galaxy. 


Rotation measure search from FAST observations A search for rotation measure was per- 
formed with the FAST data. The polarization was calibrated by correcting for differential gains 
and phases between the receptors through separate measurements of a noise diode injected at an 
angle of 45° from the linear receptors. We searched for the RM from —3.0x 10? to 3.0x 10? rad m ?. 
No significant peak was found in the Faraday spectrum. The observed lack of polarization could 


be due to the intrinsic low linear polarization, a depolarization process within the source or from 


intra-channel Faraday rotation. 
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Extended Data Legends 


ED Figure 1 | Positions of the bursts and persistent radio source identified with VLA 
observations. The observation were performed at 1.5, 3, and 5.5 GHz. They were shown as 
offsets from the best-fit position of the ensemble of bursts, at R.A. = 16h02m04.272s, Decl. = 
—11?17'17.32"(J2000). The uncertainties on the positions of the bursts are indicated with shaded 
ellipses, and those for the PRS are shown with error bars. These uncertainties include 1c statistical 


errors and estimates for systematic errors added in quadrature. 


ED Figure 2 | Spectrum of the PRS associated with FRB 20190520B. The bandwidth 
is split into two subbands for all the observations and the corresponding radio flux densities are 
shown. The observation dates in MJD are shown in the legend. The frequencies for each subband 
are shifted slightly in order to show the flux density error bars. The inset plot shows the average 
flux at each subband. A power-law model fit to these measurements is shown by the black dashed 
line and yields a spectral index of —0.41 + 0.04 for the PRS. 


ED Figure 3 | Optical spectrum of the FRB 20190520B host galaxy obtained at Palo- 
mar. The redshift of z = 0.241 was determined using the [OIII]-5007À line and Ha line (> 5c 
detections). These two emission lines are narrow, with FWHM - 10 A. The flux scale of the 


spectrum is not corrected for slit loss or extinction. 


ED Figure 4 | Posterior probability distributions. The probability distribution are for the 
shape parameter k and event rate r of the Weibull distribution. a) Results from a fit to all burst 
waiting times of FRB 20190520B detected by FAST. b) Results from a fit to all burst waiting times 
longer than 1 s. These bursts are at fluences higher than 9 mJy ms. 


ED Figure 5 | The fluence-width distribution at 1.25 GHz for FRB 20190520B. The 
widths and fluences of the 79 bursts detected by FAST are also shown in supplementary informa- 
tion Table 1. 


ED Figure 6 | The distribution of burst energies for FRB 20190520B. The bursts are 
from FAST observations. The red line and black dashed line are the log-normal fit and 9096 
completeness threshold, respectively. The 9096 completeness threshold uses 0.023 Jy ms, which is 
the simulation result from ref. 24 


ED Figure7 | Dynamic spectra of bursts with and without scattering. a) Frequency- 
time dynamic spectra of bursts with significant evidence of pulse broadening, along with 1D burst 
profiles that have been averaged across the entire frequency band. b) Examples of bursts without 


significant evidence of pulse broadening. White lines indicate excised radio frequency interfer- 
ence. All 1D burst profiles are shown in units of the signal-to-noise ratio (S/N). For P31, there is 
a potential combination of scattering, time-frequency drift, and multiple unresolved burst compo- 
nents; in this case we report a scattering time that should be interpreted as an upper limit. Bursts 


with scattering time constraints are shown in supplementary information Table 1. 


ED Table 1 | VLA observations of the PRS and number of bursts detected in each 
observation. Frequency refers to the center of the observing band and duration represents the total 
length of on-source observations. Note that the realfast system was not run on the last observation 


due to a system error. 


ED Table 2 | Localised positions of the PRS from 1.5, 3, 5.5 GHz VLA deep images. 
The first column shows the observing central frequency. The second and third columns report the 
coordinates of the PRS from the deep images. The remaining columns show the 1 c statistical 


errors and the cross-matched systematic offsets in right ascension and declination. 


ED Table 3 | Localised positions of the VLA bursts. The label in the first column shows 
the frequency band of observation, followed by the burst number. DM represents the S/N maxi- 
mizing DM obtained using offline refinement of the bursts and S/N reports the maximum image 
S/N. The remaining columns show the burst positions (and errors) obtained using CASA calibra- 
tion and image fitting. The systematic errors reported here were estimated using deep radio images 


generated using all the observations at the respective frequency band. 


ED Table 4 | Properties of the VLA bursts. The label in the first column shows the fre- 
quency band of observation, followed by the burst number. MJD is referenced to the solar system 
barycenter (in the TDB scale) and corrected for dispersion delay at infinite frequency. Width values 
here should be considered as upper limits. Frequency represents the range of frequencies in which 
the burst signal is prominently seen, and was visually determined. The DM that maximizes the S/N 
of the respective burst is given in the third column. The apparent DMs may reflect time-frequency 
structure, as opposed to bona fide DM variations. 


Extended Data Legends 


0.6 
0.4 
0.2 
~ 0.0 
o 
Oo 
8 
a 
È -0.2 
< 
-0.4 
-0.6 
e 1.5 GHz Bursts 
* 3 GHz Bursts 
e 5.5 GHz Bursts 
-0.8 Æ PRS (1.5 GHz) 
Æ PRS (3 GHz) 
Æ PRS (5.5 GHz) 


-0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 
A RA (arcsec) 


Extended Data | Figure 1 : Positions of the bursts and persistent radio source identified with VLA ob- 
servations. The observation were performed at 1.5, 3, and 5.5 GHz. They were shown as offsets from the best-fit 
position of the ensemble of bursts, at R.A. = 16h02m04.272s, Decl. = —11?17'17.32"(J2000). The uncertainties on 
the positions of the bursts are indicated with shaded ellipses, and those for the PRS are shown with error bars. These 
uncertainties include lo statistical errors and estimates for systematic errors added in quadrature. 
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Extended Data | Figure 2 : Spectrum of the PRS associated with FRB 20190520B. The bandwidth is split 
into two subbands for all the observations and the corresponding radio flux densities are shown. The observation dates 
in MJD are shown in the legend. The frequencies for each subband are shifted slightly in order to show the flux density 
error bars. The inset plot shows the average flux at each subband. A power-law model fit to these measurements is 
shown by the black dashed line and yields a spectral index of —0.41 + 0.04 for the PRS. 
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Extended Data | Figure 3 : Optical spectrum of the FRB 20190520B host galaxy obtained at Palomar. The 
redshift of z 2 0.241 was determined using the [OIII]-5007À line and Ha line (> 5c detections). These two emission 


lines are narrow, with FWHM - 10 À. The flux scale of the spectrum is not corrected for slit loss or extinction. 
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Extended Data | Figure 4 : Posterior probability distributions. The probability distribution are for the 
shape parameter k and event rate r of the Weibull distribution. a) Results from a fit to all burst waiting times of FRB 
20190520B detected by FAST. b) Results from a fit to all burst waiting times longer than 1 s. These bursts are at 
fluences higher than 9 mJy ms. 
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Extended Data | Figure 5 : The fluence-width distribution at 1.25 GHz for FRB 20190520B. The widths 
and fluences of the 79 bursts detected by FAST are also shown in supplementary information Table 1. 


80 


c f i ES an | 
3 60 + | 
o | ! 1 
s ^| MEE 
C L ! J 
E 20} —-— Burst of FRB 190520B | 
a | — — LogNormal fit | 
15 — — | S 
l 
1 
90% completeness | 
zol threshold l Va | 
È 
T 
= \ 
fo) I 
d 
5r i J 
I 
EE. | 
10° 108° 109? 


Energy (erg) 


Extended Data | Figure 6 : The distribution of burst energies for FRB 20190520B. The bursts are from FAST 
observations. The red line and black dashed line are the log-normal fit and 9096 completeness threshold, respectively. 
The 90% completeness threshold is 0.023 Jy ms. 
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Extended Data | Figure 7 : Dynamic spectra of bursts with and without scattering. a) Frequency-time 
dynamic spectra of bursts with significant evidence of pulse broadening, along with 1D burst profiles that have been 
averaged across the entire frequency band. b) Examples of bursts without significant evidence of pulse broadening. 
White lines indicate excised radio frequency interference. All 1D burst profiles are shown in units of the signal-to- 
noise ratio (S/N). For P31, there is a potential combination of scattering, time-frequency drift, and multiple unresolved 
burst components; in this case we report a scattering time that should be interpreted as an upper limit. Bursts with 
scattering time constraints are shown in supplementary information Table 1. 


Start Time Frequency | Duration | Beam size | Flux density | Bursts 
(MJD) (GHz) (min) C.) (uy) 
59051.033815 1.5 96 4.78x2.90 25816 1 
59053.046790 1.5 96 4.63x2.90 27323 2 
59079.004955 5.5 41 1.38x1.02 145+10 0 
59079.971533 5.5 41 1.44x1.05 164411 0 
59090.941498 5.5 41 1.45x1.04 158x9 0 
59091.938768 3 41 2.43x1.85 195+14 5 
59104.867225 3 41 2.76X1.76 160+13 0 
59104.908789 5.5 41 1.43x1.04 15129 0 
59105.976837 3 41 2.44x1.60 186+15 0 
59107.915546 5.5 41 1.38x1.02 153+10 1 
59111.116215 3 41 5.00x1.45 176+16 0 
59161.677240 5.5 41 1.92x0.47 13913 0 
59167.637761 3 41 2.76x0.51 233417 0 
59169.640486 3 41 2.47x0.48 211414 - 


Extended Data | Table 1 : VLA observations of the PRS and number of bursts detected in 
each observation. Frequency refers to the center of the observing band and duration represents 
the total length of on-source observations. Note that the realfast system was not run on the last 


observation due to a system error. 


Frequency R.A. Decl. Statistical Error Systematic Offsets 
(GHz) R.A. (^) Decl. (”) R.A. (^) Decl. (”) 
1.5 16h02m04.2543s  —11dl7ml17.4146s 0.0386 0.0713 0.037340.2460 0.0144+0.3035 
3 16h02m04.2611s  —11d17m17.3869s 0.0176 0.0168  —0.06964x0.2772 0.0382+0.1552 
5:5 16h02m04.2618s  —11d17m17.3375s . 0.0075 0.00905 -—0.0947+0.1161 0.1069+0.0554 


Extended Data | Table 2 : Localised positions of the PRS from 1.5, 3, 5.5 GHz VLA deep 
images. The first column shows the observing central frequency. The second and third columns 
report the coordinates of the PRS from the deep images. The remaining columns show the 1 c 


statistical errors and the cross-matched systematic offsets in right ascension and declination. 


Band S/N R.A. Decl. Statistical Error Systematic Error 
R.A.(") Decl. (~) R.A.) Decl. () 
L1 16  16h02m04.27428 -—11d17m17.0535s 0.124 0.208 0.2460 0.3035 
L2 27  16h02m04.2847s -—11d17m17.1912s 0.060 0.145 0.2460 0.3035 
L3 15  16h02m04.2691s  —11dl7ml17.1299s 0.070 0.139 0.2460 0.3035 
S1 19 16h02m04.2567s —11d17m17.4093s 0.044 0.080 0.2772 0.1552 
S2 17 16h02m04.2449s = —11d17m17.1874s 0.153 0.173 0.2772 0.1552 
S3 15  16h02m04.2602s  —11dl7ml17.5078s 0.103 0.176 0.2772 0.1552 
S4 10 16h02m04.2661s  —11dl17ml7.3666s . 0.333 0.247 0.2772 0.1552 
S5 7 16h02m04.2500s  -11dl17m17.9337s | 0.214 0.269 0.2772 0.1552 
Cl 19  16h02m04.2734s  —11dl7ml17.3078s — 0.041 0.071 0.1161 0.0554 


Extended Data | Table 3 : Localised positions of the VLA bursts. The label in the first column 
shows the frequency band of observation, followed by the burst number. S/N reports the maximum 
image S/N. The remaining columns show the burst positions (and errors) obtained using CASA 
calibration and image fitting. The systematic errors reported here were estimated using deep radio 


images generated using all the observations at the respective frequency band. 


Name MJD DM Flux density (Jy) Width (ms) Frequency (GHz) 
L1 59051.0476715 1209.60 0.113+0.006 20 1.33-1.52 
L2 59053.1045803 1227.15 0.275+0.010 20 1.39-1.52 
L3 59053.0973272 1236.60 0.189+0.009 10 1.78-1.97 
S1 59091.9604318 1291.30 0.075+0.004 20 3.12-3.50 
S2 59091.9560792 1193.40 0.100+0.009 20 2.63-2.88 
S3 59091.9635008 1222.75 0.093+0.007 20 2.50-2.76 
S4 59091.9435287 1216.80 0.037+0.005 10 2.88-3.24 
S5 59091.9668476 1276.50 0.022+0.003 10 2.50-3.50 
C1 59107.9272138 1267.50 0.068+0.004 10 4.49—4.88 


Extended Data | Table 4 : Properties of the VLA bursts. The label in the first column shows the 
frequency band of observation, followed by the burst number. MJD is referenced to the solar sys- 
tem barycenter (in the TDB scale) and corrected for dispersion delay at infinite frequency. Width 
values here should be considered as upper limits. Frequency represents the range of frequencies 
in which the burst signal is prominently seen, and was visually determined. The DM that maxi- 
mizes the S/N of the respective burst is given in the third column. The apparent DMs may reflect 


time-frequency structure, as opposed to bona fide DM variations. 


Supplementary Information | Table 1 : Properties of the 79 FRB 20190520B bursts detected by FAST. 


Burst ID^ Burst time” DM* Pulse width Scatter tail? Fluence Energy 
(MJD) ( pe cm~?) (ms) (ms) (mJy ms) (x10 erg) 

P1 58623.716958197 1238.6(20) 8.7(8) - 75(3) 9.1(4) 
P2 58623.716997034 . 1238.6(20) 16.4(6) - 30(2) 3.6(2) 
P3 58623.717075743 . 1238.6(20) 7.1(8) - 53(3) 6.4(4) 
P4 58623.717421629  1238.6(20) 7.3(10) - 51(4) 6.2(5) 
P5 58963.766862939  1209.1(10) 10.4(2) - 65(6) 8.0(8) 
P6 58963.790991383 . 1209.1(10) 10.2(17) 13.123) 52(5) 6.4(6) 
P7 58991.684351079 . 1214.2(5) 20.3(4) - 107(18) 13.022) 
P8 58991.686137870 . 1214.2(5) 22.8(83) - 100(10) 12.2(12) 
P9 58991.686875951 1214.2(5) 9.2(2) - 81(12) 9.8(15) 
P10 58991.704636888 1214.2(5) 19.5(17) 3.6(8) 108(12) 13.1(15) 
P11 58991.704638485 1214.2(5) 9.7(7) 42(0) 5.2(0) 


P12 58991.717697335 . 1214.2(5) 11.9(11) 10.3(12) 12221) 14.9(25) 
P13 58991.717878516 — 1214.2(5) 14.3(40) 5.1(25) 58(10) 7.1(12) 
P14 58991.718217265 . 1214.2(5) 11.8(15) - 63(6) 7.I(8) 

P15 58991.718217670 . 1214.2(5) 17.5(3) - 39(4) 4.8(5) 

P16 58991.718218082 . 1214.2(5) 20.2(22) 12.6(19) 153(28) 18.6(34) 
P17 58991.735450006 . 121422(5) 30.2(5) - 138(21) 16.8(26) 
P18 58991.750640038 . 121425) 17.7(13) - 60(9) 7.4(11) 
P19 58991.750640790 = 1214.2(5) 20.2(5) - 81(14) 9.9(17) 
P20 59060.484475154 — 1209.1(5) 13.5(13) 11.3(14) 241(12) 29.3(15) 
P21 59060.507858658 . 1209.1(5) 7.3(9) 9.3(12) 170(8) 20.7(10) 
P22 59060.525960600 . 1209.1(5) 17.3(24) 7.6(16) 167(14) 20.3(17) 
P23 59061.512755570 — 1209.1(6) 8.5) - 486) 5.907) 

P24 59061.512755780 . 1209.1(6) 10.0(23) 4.7(17) 129(16) 15.7(20) 


P25 59061.516277966  1209.1(6) 7.1(6) 29(2) 3.6(3) 
P26 59061.516278904 1209.1(6) 7.5(4) : 49(5) 6.0(7) 
P27 59061.516279436 1209.1(6) 11.4(8) : 49(4) 6.0(5) 
P28 59061.524341261 1209.1(6) 10.7(8) 10.4(9) 122(16) — 14.8(19) 


P29 59061.535633328  1209.1(6) 14.0(4) 84(8) 10.2(10) 
P30 59061.536568858  1209.1(60 21.907) . 72(9) 871) 
P31 59061.536569628 — 1209.1(6) 19.0014) 14.03) — 122(4)  14.8(50) 
P32 59061.537903820  1209.1(600  23.9(27) 4.9(13) 65(8) 8.0(9) 

P33 59061.539298603 1209.1(6) 9.9(3) : 73(9) 8.9(11) 
P34 59061.541828100  1209.1(6)  19.8(20) 87(10) — 2128) 25.834) 
P35 59067467544345 1202.86) 15.2(76) : 67(8) 8.2(10) 
P36 59067.486738799 1202.8(6)  14.9(22) 17.829 — 131(14) 15917) 
P37 559067.486739378 1202.8(6) 33.1(7) à 146(18) — 17.821) 
P38 59067.502691880  1202.8(6) 14.3(6) : 134(3) 16.3(3) 
P39 59067.500899127  1202.8(60 12218) 10.3(17) 122(8) 14.8(9) 
P40 59067.535246460  1202.8(6) 7.62) : 74(9) 9.0(11) 
P41 59069.495909561  1190.2(11) 16.865) : 14300) 17402) 
P42 59069.501196109 1190211) 18.7(1) a 291(14) — 354(17) 
P43 59060.514994796  11902(1))  12.2(15) 11.617) 192(13) 23415) 


P44 59071.472522775 1200.001) 52(3) : 79(1) 9.6(2) 
P45 . 59071.472523007  1200.0(11) 7.2(7) , 66(1) 8.1(2) 

P46 59071.491696655 1200.01) 10.4(20) 4.9(15) 184(4) 22.4(5) 
P47 59073496887082 1197007) 7.8Q) z 10709) BIOD 
P48 59073.515256071  1197.0(7) 18.7(6) 3 147(2) 17.9(2) 
P49 59073.515256881 1197.0(7) 11.1(7) - 57(4) 7.0(5) 

P50 59075.454353002  1210.6(6) 249(1) z 31319) 38003) 
P51 59075.454862469 — 1210.6(6) 22.0(1) : 20912) — 25.5(15) 
P52 59075472181012  1210.6(6) 14.0(5) : 114(6) 13.8(8) 
P53 59075.484186304 1210.6(6) 6.0(0) : 15706) — 19.120) 
P54 59075.496463775 1210.6(6) 16.7(4) x 83(8) 10.19) 


P55 59077.448938437 . 1209:9(4) IT.9(5) - 103(7) 12.5(8) 


Supplementary Information | Table 1 : Properties of the 79 FRB 20190520B bursts detected by FAST. 


Burst ID^ Burst time” DM* Pulse width Scatter tail’ Fluence Energy* 
(MJD) ( pc cm~?) (ms) (ms) (mJy ms) (x10 erg) 
P56 59077.448939131 1209.9(4) 14.4(7) - 107(5) 13.1 


P57 59077.449538432 1209.9(4) 13.4(17) 10.6(16) 199(10) 24.2(12) 
P58 59077.449939184 . 1209.9(4) 5.4(26) - 148(7) 18.0(8) 
P59 59077.449939372  1209.9(4) 9.6(29) - 121(11) 14.7(14) 
P60 59077.460100338 . 1209.9(4) 13.0(23) 13.5(28) 132(1) 16.1(2) 
P61 59077.460503968 . 1209.9(4) 10.3(18) 10.6(21) 94(5) 11.4(6) 
P62 59077.466295203 . 1209.9(4) 11.7(2) - 156(10) 19.0(12) 
P63 59077.468715062  1209.9(4) 18.4(35) 7.6(18) 128(9) 15.6(11) 
P64 59077.469903193 . 1209.9(4) 11.2(10) 11.3(12) 152(7) 18.5(9) 
P65 59077.475047947 . 1209.9(4) 11.0(17) 8.3(16) 107(6) 13.0(8) 


P66 —-559077.475330491  1209.9(4) 14.0(8) 12.99) — 205101)  249(14) 
P67  59077.477447949  1209.9(4) 2.0(1) : 70(2) 8.5(3) 
P68 59077485451841  1209.9(4) 14.3(1) : 332(7) 40.4(8) 
P69 59077485451957  1209.9(4) 14.4(1) : 262005) — 31.8(19) 
P70 59077.490027004  1209.9(4) 8.4(4) : 132(2) 16.1(3) 


P71 59077.491413589 . 1209.9(4) 10.3(3) - 115(5) 14.0(6) 
P72 59077.497806611 1209.9(4) 15.3(22) 10.7(21) 138(9) 16.8(11) 
P73 59077.497957805 = 1209.9(4) 13.6(4) 89(5) 10.9(6) 


P74 59077.4989609066  1209.9(4) 4.6(7) 11.8(15) 94(6) 11.4(7) 
P75 59089.428163703 1186.725) 1120) 2 112) 13.6(4) 
P76  59089.436165435  1186.7(25) 13.3(6) z 83(2) 10.1(3) 
P77 59111.370525950 1183.307) 20.6(6) z 9911) 12.0(14) 
P78 . 59111.370526098  1183.3(17) 6.7(1) : 106(8) 12.9(9) 
P79 59111.370978062  1183.3(17) 14.6(5) . 98(9) 11.9(11) 


“ The burst IDs are from 1 to 79. P1 to P4 were detected in the drift 
scan mode. P4 to P79 were detected in the tracking mode. 

^ Arrival time of burst at the solar system barycenter in barycentric 
dynamical time (TDB). They are corrected to the frequency of 1.5 
GHz in the International Celestial Reference System (ICRS). 

* All the bursts detected on the same day are assigned the best fit DM 
value of the highest S/N burst from that day; the highest S/N burst ID 
is shown in bold for each observing epoch. Epochs are separated by 
single horizontal lines and the apparent DM variations may solely be 
the result of the variable spectra-temporal structure of the bursts. 

d We only report burst scattering timescales with fractional uncertain- 
ties less than 5096. These scattering times are based on fitting a Gaus- 
sian pulse convolved with a one-sided exponential to the 1D burst pro- 
file, which was averaged over the entire burst bandwidth. 

* Energy here refers to equivalent isotropic energy. 

f Reported uncertainties correspond to the last consecutive digits; e.g., 
1238.6(20) = 1238.6 + 2.0. 


